Ultrahigh energy neutrinos and non-hnear QCD dynamics 



Magno V.T. Machado '^•''^^ 
°- High Energy Physics Phenomenology Group, GFPAE IF-UFRGS 
GP 15051, GEP 91501-970, Porto Alegre/RS, Brazil 
' The Ahdus Salam International Centre for Theoretical Physics (ICTP) 
Strada Costiera 11, 34014 Trieste, Italy 
" GERN TH Division, CH-1211 Geneve 23, Switzerland 

The ultrahigh energy neutrino-nucleon cross sections are computed taking into account different 
phenomenological implementations of the non-linear QCD dynamics. Based on the color dipole 
framework, the results for the saturation model supplemented by DGLAP evolution as well as for 
the BFKL formalism in the geometric scaling regime are presented. They are contrasted with recent 
calculations using NLO DGLAP and unified BFKL-DGLAP formalisms. 
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I. INTRODUCTION 

The inelastic interaction of neutrinos with nucleons is 
in general described by the QCD-improved parton model. 
In ultrahigh energy neutrino interactions one probes a 
kinematical region which is not accessible in the current 
collider experiments. It is known so far that these reac- 
tions can be sensitive to the domain for the usual Bjorken 
variable x ~ m^^^/S^ ^ lO'^ at - lO^^ GeV 
and electroweak boson virtualities ~ ""^w z ~ 
GeV^, where mw,z are the boson masses. In contrast, 
in the current electron-proton collider DESY HERA the 
interactions are probed in the typical kinematical range 
X > 10-5 and < < 10^ GeV^ (the variables x 
and photon virtuality are correlated, with small x 
data associated to low virtualities). There, the mea- 
surements of the deep inelastic ep scattering structure 
functions have put important constraints on parton dis- 
tributions, mainly sea-quarks and gluons, in the small x 
region. In perturbative QCD these distributions are ex- 
pected to grow as x decreases, which has been strongly 
confirmed by HERA data. Therefore, in order to provide 
accurate predictions for the ultrahigh energy neutrino- 
nucleon cross sections, a precise extrapolation of the 
structure functions to the small x and large region 
probed in these interactions is needed. 

In the usual QCD-improved parton model, the power 
growth of the parton distributions functions in very high 
neutrino energies lead consequently to a power increasing 
of the neutrino-nucleon total cross section, which in turn 
implies in a violation of unitarity at high energies. This 
issue is an outstanding theoretical challenge, which has 
produced intense work towards a complete understand- 
ing of the non-linear (higher twist) QCD corrections in 
the high parton density regime at fixed impact parame- 
ter (for a recent review, see [l| and references therein). 
At present, one has reliable estimates for the kinemati- 
cal domain where non-linear (saturation) effects should 
play an important role. Namely, below a typical trans- 
verse momentum scale Qsat(^) x~^, the so-called sat- 
uration scale the growth of parton (gluon) density 



is tamed towards the black-disk limit of the target. This 
slow down of the gluon distribution function prevents the 
total cross section violating unitarity. An important pre- 
diction of the non-linear QCD approaches is the property 
of geometric scaling jj] below saturation scale, where the 
gluon density (and cross sections) scales with the satura- 
tion momentum, which grows with l/a: as a power. It has 
also been shown that such scaling survives even above 
the saturation regime ^ IS- the kinematical region 
has been estimated to be Ql^^{x^ < < Q'^^^{x) / Aq^^^ 
relying on the BFKL equation 7] in the presence of sat- 
uration This feature has important consequences in 
the underlying QCD dynamics at the nucleus-nucleus ac- 
celerators RHIC and LHC and should pose strong con- 
straints to the cross sections beyond current accelerator 
energies. Concerning high energy neutrino cross sections, 
there are estimates constrained by HERA data based 
on DGLAP M and effects of BFKL hnear evolutions 
0, |lO|, [llj, I12II13I Hij as well as implementations of non- 
linear QCD corrections The linear NLO DGLAP 
and unified BFKL-DGLAP results are found to be con- 
sistent with each other presenting somewhat small devi- 
ations 01 . An interesting analytical calculation using 
an approximate DGLAP solution with initial conditions 
from a soft non-perturbative model has been recently re- 
ported in Ref. 16]. Saturation effects are shown to be 
quite small in recent DGLAP analysis 0|, whereas it 
was found a reduction in the cross section by a factor 2 
at El, ~ 10^^ GeV considering the saturation model or 
screening effects via Balitsky-Kovchegov equation [Tsj . 
Recently, it has been claimed in Ref. ^3 that the ge- 
ometric scaling property can lead to an enhancement of 
the neutrino-nucleon total cross section by an order of 
magnitude compared to the leading twist cross section 
and would lead it further to its unitarization. 

In this work one presents estimates for the ultrahigh 
energy neutrino-nucleon cross sections, focusing on the 
saturation models which encode the physics of non-linear 
QCD evolution in a phenomenological way. The main 
point is those models provide a suitable dynamical ex- 
trapolation to very high energies/virtualities in very con- 
trast with usual DGLAP based approaches, which are 
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not a priori valid in that regime and are dependent on 
quark and gluon distributions also unknown (no con- 
straint from accelerator data) in that kinematical limit. 
Based on the color dipole framework, first one presents 
the result for the saturation model 18] supplemente d by 
QCD evolution The original saturation model [isj 
was proposed for either low values of < 150 GeV^ , 
whereas high energy neutrino interactions are sensitive to 
large virtualities ~ "i^wz- Therefore, it is necessary 
to perform calculations with complete QCD evolution. 
Furthermore, one computes estimates of the ultrahigh 
neutrino-nucleon cross section through an implementa- 
tion of the dipole cross section considering BFKL for- 
malism in the geometric scaling region 20]. In the last 
section, the numerical results for both cases are shown 
and contrasted with recent NLO DGLAP estimates and 
with unified BFKL-DGLAP formalism without satura- 
tion effects. 



II. NEUTRINO-NUCLEON CROSS SECTION 

Deep inelastic neutrino scattering can proceed via 
or Z° exchanges. The first case corresponds to charged 
current (CC) and the second one to the neutral current 
(NC) interactions, respectively. The standard kinemati- 
cal variables which describe the processes above are given 
by, s = 2mNE^, = 



2; = ^ and y 



2p-q ^^^^ y jtinEj^ ' 

Here, rriN is the nucleon mass, E^, labels the neutrino en- 
ergy and p and q are the four momenta of the nucleon and 
of the exchanged boson, respectively. The usual Bjorken 
variable and the momentum transfer are usually labeled 
as X and Q^, whereas s is the total centre-of-mass energy 
squared. One assumes an isoscalar target N — {p + n)/2, 
which is a good approximation for the present purpose. 

At high energies a quite successful framework to de- 
scribe QCD interactions is provided by the color dipole 
formalism |22j] . which allow an all twist computation (in 
contrast with the usual leading twist approximation) of 
the structure functions. The physical picture of the inter- 
action is the deep inelastic scattering at low x viewed as 
the result of the interaction of a color qq dipole which the 
gauge bosons fluctuate to with the nucleon target. The 
interaction is modeled via the dipole-target cross section, 
whereas the boson fluctuation in a color dipole is given 
by the corresponding wave function. The DIS structure 
functions for neutrino-nucleon scattering in the dipole 
picture read as [T^IIitI]. 
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where r denotes the transverse size of the color dipole, z 
the longitudinal momentum fraction carried by a quark 
and iPt^'l proportional to the wave functions of the 
(virtual) charged or neutral gauge bosons correspond- 
ing to their transverse or longitudinal polarization (F2 ~ 
Ft + Fl). Explicit expressions for the boson {W^ and 



Z") wave functions squared are as follows pjl ], 

^[z' + {l-zf]e^Kf{er), (2) 
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)1^ = -z^il-zYQ^K^ier), (3) 

IV'f (r,z,Q2)|2 ^ ^Kwz'{\~zfQ'Kl{eT), (5) 

where one defines the auxiliary variable e — z(\ ~ z) 
and Ko^i{x) are the Mc Donald's functions. Here, Kw = 
(L^ + + R'^ + B?j) and the the chiral couplings are ex- 
pressed as functions of the Weinberg angle 9w as follows. 
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Following Ref. here one considers only four fla- 

vors {u,d,s,c) and assumes them massless, whereas it 
has been shown Q that heavy quarks {b,t) give rela- 
tively small contribution. Moreover, the color dipoles 
contributing to Cabibbo favored transitions are ud (du) , 
cs (sc) for CC interactions and uu, dd, cc, ss for NC inter- 
actions. In our further numerical calculations using the 

CC NC 

dipole framework, we use the structure function F^ ' 
given by the usual LO DGLAP expressions [IJ], despite 
their contributions to be small for the present purpose. 

The total CC (NC) neutrino-nucleon cross sections as 
a function of the neutrino energy are given by the inte- 
gration over available phase space at the given neutrino 
energy. They read as. 
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where y = Q'^ /{xs) and a minimum value (5^i,j (of order 
of a few GeV's) on is introduced in order to stay in 
the deep inelastic region. 

In the present work we are interested in the ultrahigh 
energy neutrinos {Ei, 3> 10 GeV), where the valence 
quark contribution stays constant and physics is driven 
by sea quark contributions. In this region, we will use the 
available phenomenology on non-linear QCD evolutions 
to study saturation effects in the neutrino-nucleon cross 
section. 

The dipole cross section adip, describing the dipole- 
nucleon interaction, is substantially affected by satura- 
tion effects at dipole sizes r > l/Qsat- Here, we fol- 
low the saturation model , which interpolates between 
the small and large dipole configurations, providing color 
transparency behavior, adip , as r — > and constant 
behavior, a^ip ~ do, at large dipoles. The transition is 
rendered by the saturation phenomenon. The parametri- 
sation for the dipole cross section takes the eikonal-like 
form. 
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where the saturation scale Q^{x) = (^) defines the 
onset of the saturation effects. As known, most of the 
contribution to ultrahigh energy neutrino cross section 
comes from very small x w my^ ^ /2mjqEi, and therefore 
a sizeable part of the contributions is in the geometric 
scaling region or even in the saturation region The 
parameters (cto, A and xq) were obtained from a fit to 
the HERA data The variable i = a; ( 1 + 4m^/Q2) 

gives a suitable transition to the photoproduction re- 
gion. The saturation model has been used to compute 
neutrino- nucleon cross section in Ref. |l5j| . There, it 
was shown that the saturation model does not include 
complete DGLAP QCD evolution at high virtualities and 
tends to be therefore lower than other calculations, prob- 
ably meaning that it is rather incomplete and inaccurate 
in the high region. 

Recently, a new implementation of the model includ- 
ing QCD evolution (labeled BGBK model) has ap- 
peared. Now, the dipole cross section depends on the 
gluon distribution as, 
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where the initial condition at /i^ = 1 GeV^ is x G = 
Ag x~'^<' (1 — x)^-^ and = C/r^ + ^q. The phenomeno- 
logical parameters are determined from a fit to small x 
HERA data. The function G{x,fi^) is evolved with the 
leading order DGLAP evolution equation for the gluon 
density with initial scale Qq — 1 GeV^ . The improvement 
preserves the main features of the low-Q^ and transition 
regions, while providing QCD evolution in the large-Q^ 
domain. 

Despite the saturation model to be very successful in 
describing HERA data, its functional form is only an ap- 
proximation motivated by the Glauber-Mueller formula, 
which it does not include impact parameter dependence. 
On the other hand, an analytical expression for the dipole 
cross section can be obtained within the BFKL formal- 
ism. Currently, intense theoretical studies has been per- 
formed towards an understanding of the BFKL approach 
in the border of the saturation region ^ |^ . In particu- 
lar, the dipole cross section has been calculated in both 
LO and NLO BFKL ^ approach in the geometric 
scaling region 24]. It reads as. 



cTdip = Cto [r^Q^^^ix)] exp 
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where cto = 2ttR^ {Rp is the proton radius) is the overall 
normalization and the power 7sat is the (BFKL) saddle 
point in the vicinity of the saturation line — Qgg^^{x) 
(the anomalous dimension is defined as 7 = 1 — 7sat). 
As usual in the BFKL formalism, ctg = Ncas/ir, (3 ~ 
28C(3) and the notation Y — ln(l/a;). The quadratic 
diffusion factor in the exponential gives rise to the scaling 
violations. Eq. has been used in Ref. ^3 to show 
an estimation of the enhancement factor (about a factor 



10) in ultrahigh neutrino- nucleon cross section in contrast 
with the usual leading twist calculations. 

The dipole cross section in Eq. does not in- 

clude an extrapolation from the geometric scaling re- 
gion to the saturation region. This has been recently 
implemented in Ref. |2n| |. where the dipole amplitude 
A/'(x, r) = adip/2TTRp was constructed to smoothly inter- 
pole between the limiting behaviors analytically under 
control: the solution of the BFKL equation for small 
dipole sizes, r <^ 1/Qsat(a^), and the Levin- Tuchin law 
psf for larger ones, r l/Qsatix). A fit to the struc- 
ture function F2(x, Q^) was performed in the kinematical 
range of interest, showing that it is not very sensitive to 
the details of the interpolation. The dipole cross section 
was parametrized as follows. 
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where the expression for rQsatix) > 2 (saturation region) 
has the correct functional form, as obtained either by 
solving the Balitsky-Kovchegov (BK) equation [2^, or 
from the theory of the Color Glass Condensate (CGC) 
Hereafter, we label the model above by CGC. The 
coefficients a and b are determined from the continuity 
conditions of the dipole cross section at rQ sat{x) — 2. 
The coefficients 7sat = 0.63 and k = 9.9 are fixed from 
their LO BFKL values. In our further calculations it 
will be used the parameters Rp = 0.641 fm, A = 0.253, 
xo = 0.267 X 10-"^ and Mo = 0.7, which give the best fit 
result. A large x threshold factor (1 — x)^ will be also 
considered. 



III. RESULTS AND DISCUSSION 

In what follows, one computes the ultrahigh energy CC 
and NC neutrino-nucleon cross sections using the differ- 
ent models containing saturation effects. In Fig. ^one 
presents the results for the saturation model (GBW), the 
BGBK model and also one includes a DGLAP calcula- 
tion. The result for the saturation model (dot-dashed 
lines) has been recently computed in Ref. (15| , which we 
reproduce here. The introduction of the DGLAP evolu- 
tion in the saturation model enhances the cross section by 
a factor 3 or even more at E,y ~ 10"'^^ GeV. This is shown 
by the BGBK curves (solid lines). It should be noticed 
that the power growth for saturation model is milder 
than the BGBK model and it gives an upper limit for 
the role played by saturation effects in the cross section. 
It should be emphasized this assertion concerns to a com- 
parison between phenomcnological saturation model and 
its version including DGLAP evolution. Namely, other 
non-linear QCD approaches can give a larger amount of 
saturation in the same kinenatical regime and the phe- 
nomcnological saturation model probably is inaccurate 
in this region. For sake of completeness, the DGLAP re- 
sults (LO and NLO) are also presented. For this purpose. 
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FIG. 1: The v N charged current (a) and neutral current (b) 
as a function of neutrino energy. The curves correspond to 
the results for the GBW saturation model (dot-dashed curves) 
and the BGBK model (solid lines). The LO and NLO DGLAP 
results are represented by dashed and long-dashed curves, re- 
spectively. 



FIG. 2: The v N charged current (a) and neutral current (b) 
as a function of neutrino energy. The solid curves correspond 
to the results for the CGC model. The NLO DGLAP resuH 
is represented by the dashed and the unified BFKL-DGLAP 
formalism by the dot-dashed curves, respectively. 



we follow the calculations in Ref. where a detailed 
analysis of the NLO corrections to neutrino- nucleon cross 
sections was performed. The BGBK model has an energy 
dependence at ultrahigh energy neutrinos similar to the 
NLO DGLAP (long-dashed Hues) for > 10^ GeV. 

In order to investigate the BFKL physics in the geo- 
metric scaling region, we use the recent fit to the dipole 
cross section presented in Ref. 20J (the CGC model), 
which encodes the BFKL formalism in the geometric scal- 
ing regime and an interpolation to the saturation domain, 
as presented in the last section. It should be noticed 
that an analysis including DGLAP evolution in the fit 
would be timely. However, as the CGC model produces 
an "effective" anomalous dimension, i.e. jeS {x,rQsa,t) = 
Isat + (KAr)-Mn(4/r2Q2^J far from the geomet- 
ric scahng region r'^Ql^^{x) <^ 1, which is closer to the 
DGLAP one, the estimate presented here would be not 
strongly spoiled by QCD evolution. The results for the 
CGC model (solid lines) is presented in Fig. |21 con- 
trasted with the NLO DGLAP calculation (dashed lines) 
and also with the unified BFKL-DGLAP formalism (dot- 
dashed lines), which embodies non-leading ln(l/a;) con- 
tributions jlJI- The latter incorporates both the ln(l/a:) 
BFKL resummation and the complete LO DGLAP evo- 
lution, with the dominant non-leading ln(l/x) contri- 
butions resummed to all orders. The result found for 
CGC model is quite similar to both DGLAP and uni- 
fied BFKL-DGLAP formalisms in the range of energy 
considered here. It is not observed any enhancement in 
the neutrino-nucleon cross section in comparison with the 
leading twist DGLAP calculation, even at very high en- 
ergies of order Ey ~ 10^^ GeV. This fact is in disagree- 
ment with Ref. where it has been shown that the 
geometric scaling lead to an enhancement of the neutrino- 
nucleon total cross section by an order of magnitude ver- 



sus the leading twist cross section. We also verify that 
at El, > 10"'^^ GeV the power on energy starts to slow 
down and it should be increasingly smoother (logarith- 
mic growth) at higher energies. This expected change of 
behavior on energy of the cross section is estimated to be 
reached at E^, ^ 10^^ GeV in Ref. 0, suppressing the 
cross section compared to the leading twist result. 

We can qualitatively understand the result above an- 
alyzing the ratio between the dipole cross section for 
BFKL and DGLAP approaches in the geometric scaling 
regime. As discussed before, the ultrahigh energy neu- 
trino cross section is sensitive to the kinematical region of 



small X 



''W,Z 



/E^ and high 



Concerning 



DGLAP formalism, in such a region we are in the double 
logarithmic approximation (DLA) limit. Concerning the 
non-linear QCD approaches, at the vicinity of saturation, 
the DLA saddle point is no longe r eq ual to one, but it 
was found to be 7,^^^* ~ 1/2 HEllll. The ratio of the 
cross section in the geometric scaling region is given by. 



?sat] 
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<> sat 
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where we have disregarded the diffusion factor in the 
BFKL expression and consider ~ ~ 
The result i?scai — 1 leads to similar results between the 
all twist and leading twist calculation in the geometric 
scaling region. In fact, there is room for a slight suppres- 



sion since (7. 



BFKL 
sat 



DLA 

/sat 



0.13 and QlJrr 
On the other hand, the argument in Ref. [l3| is exactly 
true in the comparison between CGC and GBW models. 
In this case, as we will discuss in what follows, 7^f ^ = 1 
in the geometric scaling region and a sizeable enhance- 
ment between CGC and GBW is really observed. 

The qualitative analysis presented above can also help 
us to understand the enhancement of the BGBK model 



< 1. 
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(with DGLAP evolution) in relation to GBW saturation 
model. We can think in an effective "anomalous dimen- 
sion" for the saturation model in the geometric scaling 
regime: from Eq. cr^^^ ~ r'^Q'^^^^/A and then one 
verifies that 7^f^ = 1- In this case, the ratio now reads 
as. 



-Rscal 



ag.p (BGBK) 
ad,p (GBW) 



*5sat 



1 (-yf. 
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whereas at lower energies (beyond the geometric scal- 
ing regime) the DGLAP anomalous dimension takes the 
usual value and the ratio is closer to one, as expected. 

As a summary, we have presented the results for the ul- 
trahigh energy neutrino cross section considering the cur- 
rent phenomenology on non-linear QCD dynamics, which 
play an important role in the correct extrapolation to the 
regime of very high gluon density. The saturation model 
with QCD evolution was studied as well as the physics of 
BFKL dynamics at geometric scaling region. The latter 
presents suppression in the cross section at i?^ <; 10^^ 
GeV, whereas the BGBK model is consistent with NLO 
DGLAP calculations. An outstanding advantage in the 
color dipole framework is the possibility to provide ana- 
lytical expressions for the extrapolation to ultrahigh en- 



ergies, towards the black limit of the nucleon target. This 
cannot be taken into account in a leading twist approach 
as the QCD-improved parton model. Then these satura- 
tion (all twist) approaches play an important role, since 
a precise knowledge about the neutrino interactions and 
production rates is essential for estimating background, 
expected fluxes and detection probabilities [l^ l29l | . In 
particular, reliable estimates for UHE neutrino cross sec- 
tion are strongly needed for the planned experiments [s^ 
to detect them via nearly horizontal air showers in the 
Earth atmosphere, as discussed in [l^. 
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